Resonance shift effects in apertureless scanning near-field optical microscopy 



<N : 
o . 
o ■ 
(N ; 

3 ■ 
<■ 
\0 ■ 



oo 

00 

o 
o 



i 

c 

o 
o 



J. A. Porto, 1 P. Johansson, 2 S. P. Apell, 1 and T. Lopez-Rios 3 
1 Department of Applied Physics, Chalmers University of Technology and Goteborg University, S-41296 Goteborg, Sweden 
2 Department of Natural Sciences, University of Orebro, S-701 82 Orebro, Sweden 
3 Laboratoire d' Etudes des Proprietes Electroniques des Solides, 
(LEPES/CNRS), BP 166, 38042 Grenoble Cedex 9, France 
(Dated: February 1, 2008) 

We develop a theory to study apertureless scanning near-field optical microscopy which takes into 
account retardation, higher multipoles of the tip, and the multiple scattering between the tip and 
the surface. We focus on metallic systems and discuss the implication of the formation of tip-induced 
surface plasmon modes in the tip-surface system. We discuss the effects associated with the shift 
in energy of those modes as a function of the tip-surface distance. Both the local field and the 
scattering cross section are enhanced when the tip approaches the surface, but there is no general 
correspondence between the two enhancements. 

PACS numbers: 68.37.Uv,73.22.Lp,78.68.+m 



I. INTRODUCTION 

Scanning near- field optical microscopy (SNOM)lillat is 
a technique that allows to obtain optical images of ob- 
jects with suhwavelength resolution. In the pioneering 
experiments^! 2 ] subwavelength optical apertures where 
used to illuminate the sample and obtain subwavelength 
images. Nevertheless, it is also possible to perform near- 
field optical imaging without-a^tture probes. This is the 
case of apertureless SNOMjfflffijtl which has attracted 
considerable attention in recent years, allowing inter- 
esting experiments such as near-field fluorescence ijaag- 
ingja and near-field vibrational absorption detection.EHI In 
apertureless SNOM, light is focused at the tip of a scan- 
ning probe microscope (SPM) and the field enhancement 
near the tip is used to obtain the subwavelength imaging. 
Apertureless SNOM allows to perform simultaneously 
near-field optical imaging and other types of local probe 
microscopies, such as atomic force microscopy (AFM) or 
scanning tunneling microscopy (STM). In addition, it has 
provided higher resolution than other SNOM techniques. 
A key issue in the apertureless SNOM technique is the 
tip-sample coupling, and how the electromagnetic reso- 
nances at the cavity formed by the tip and the surface 
may affect the corresponding scattered field. 

The re ha s been an extensive theoretical research about 
SNOM,tilEj mostly related to imaging with aperture 
probes. With regard to apertureless SNOM, it has been 
studied theoretically considering different models j-and. 
approximations such as electric,-dipole for the tipjij'til 
the electrostatic-.approximationJl3 or the passive probe 
approximation .Ej It has also been studied in two di- 
mensionsO There have been detailed numerical calcu- 
lations r-Qf— itfhe near field associated with apertureless 
probesJl3'tj'E2l and the field enhancement haa_been dis- 
cussed in the context of aper tadess SNOM.EI In addi- 
tion, the role-,of tip geometryESEPJEJ the influence of tip 
modulation,!^ and— the implications for magneto-optical 
near-field imagingEjEZl have been studied. Nevertheless, 
there is a need for clarification of how the tip-sample 



coupling affects the signal detected and the local-field 
enhancement. In particular, the question concerning the 
link between electromagnetic resonances at the cavity 
formed by the tip and the surface and the corresponding 
scattered field needs to be addressed. These resonances 
concern free electron metals for which surface plasmons 
may exist. In the simple case of an isolated particle it is 
well known that the scattering and extinction cross sec- 
tions have resonance maxima at about the same photon 
energy. This is obvious from the fact that large charges 
at the surface give rise to large dipoles and consequently 
large radiation. If we instead consider two interacting 
particles the situation is not so simple. In this case there 
is not a direcl-jcprrelation between near-field and far-field 
resonances .E§t3 

This problem is of the highest relevance concerning 
the emerging field of optical spectroscopy with an aper- 
tureless SNOM, which comprises areas such as near-field 
Raman scattering and luminescence. The small number 
of molecules between the tip and the surface leads to ex- 
tremely small signals that are very difficult to measure. 
One possible way to overcome this problem is to take 
advantage of the resonance that may exist at the cavity 
between the tip and the surface and the extremely high 
fields associated.EII In this context, a key question is what 
is the signature of this type of resonances on the far field 
and on the optical images. 

In this article we develop a theory which allows to per- 
form a detailed study of the tip-sample interaction in 
apertureless SNOM. We model the tip by a sphere of fi- 
nite size, whose radius may correspond roughly to the 
radius of curvature of the real tip. The main features of 
the theory are: (i) retardation, (ii) higher multipoles of 
the sphere, and (iii) the multiple scattering between the 
sphere and the surface are taken into account £3 We use 
this theory to analyze the effects associated with the shift 
in energy of the surface-plasmon resonances of the tip- 
sample system as the tip approaches the surface. Those 
resonances produce an enhancement in the near field and 
in the far field as well, although there is not always an ex- 
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act correspondence between the two enhancements. We 
will also discuss the convergence of the results with the 
number of multipoles taken into account in the calcula- 
tion. As a general result, a high number of multipoles 
should be retained to arrive at convergence. 

The paper is organized as follows. In Sec. II we high- 
light the main points of the theoretical framework used 
for the calculations. In Sec. [II we discuss the numeri- 
cal results obtained for the near-field enhancement and 
scattering cross sections as a function of different param- 
eters of the system. In Sec. |v|, we present the concluding 
remarks. 
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II. METHOD 



In this section we present the. main features of the 
theoretical formalism employed.c3 Figure |l|(a) shows the 
system under study. A flat metallic sample, filling the 
half-space z < zq = —(R + d), is probed by a metallic 
tip, modeled by a sphere with radius R. The distance 
between the sphere and the surface is d, as shown in 
Fig. 0(a). 

We are interested in studying how this system scatters 
the incident light and which local electromagnetic fields 
are built up in the process. For that purpose, we consider 
a plane wave incident on this system with wave vector q, 
which can be written as (an e~ lu}t factor is omitted in 
fields and other quantities in what follows) 



E, 



t (r) = { 



(1) 



For a wave vector q, the polarization vectors for s and 
p polarization are, respectively, s = z x qii , and p = 

q x (z x qii), where is the projection of q in the 
surface plane. In the present work we will mainly focus 
on the p-polarization. 

In order to tackle the scattering off the sphere-surface 
system, the electromagnetic field inside and just outside 
the sphere is expanded in_±erms of the electric (E) and 
magnetic (M) multipoles .EJ The electric field inside the 
sphere can be expressed as 
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The electric field just outside the sphere can be written 
as the linear combination of incident and outgoing waves 
from the sphere 
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+'Vx{ [a^Ji (*t) + h t (kr)] X Im } , (3) 

where X; m is the vector spherical harmonics, ji denotes 
the spherical Bessel functions related to incoming waves, 
and hi are the spherical Hankel functions related to out- 
going waves. Equivalent expressions can be written for 
the magnetic field. 
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FIG. 1: (a) Model geometry used in the calculation and (b) 
illustration of the scheme used in the calculation 



By imposing the continuity of the tangential compo- 
nents of E and H and the normal components of B and 
D at the sphere surface, it is possible to relate the coeffi- 
cient of the outgoing wave to the incoming wave for each 
electric and magnetic multiple: 
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where s[ E ^ and s[ M ' are sphere response functions. Due 
to the symmetry of the sphere, the "quantum" number I 
is conserved in a scattering event and the response func- 
tions are independent of m. These expressions can be 
rewritten in a matrix form 



b = 



s a, 



(4) 



where *s is a diagonal tensor that relates the incoming 
wave coefficients, a, to the outgoing waves, b. 

In the situation under study, the field impinging on the 
sphere consists of three contributions [see Fig. |l](b)], 
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corresponds to the amplitude of the field from the 
original incident wave, its value can be obtained by ex- 
panding the incident plane wave in the multipoles of the 
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sphere, while a ref represents the amplitude of the incident 
wave reflected once from the sample surface. Its value is 
obtained by expanding the reflected plane wave from the 
surface in the multipolcs of the sphere. Finally, Na rep- 
resents waves scattered from the sphere that return to it 
after reflection from the sample. This term contains the 
exact field described by a, therefore it accounts for all 
multiple scattering events in which waves are scattered 
between the sphere and the sample an arbitrary number 
of times. N is a tensor in multipole space. It indicates 
how a spherical wave with a given angular momenta 1 
and m is reflected from the surface and returns to the 
sphere. In this reflection process the quantum number I 
is not necessarily conserved but, thanks to the cylindrical 
symmetry, m is. In the coupling between the sphere and 
the sample, both propagating and evanescent waves are 
taken into account. 

Retaining a finite number of multipoles, Eq. (^|) can be 
solved by a matrix inversion 
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Now, from Eq. (|j) we can obtain the coefficients for 
the outgoing waves from the sphere. From those coeffi- 
cients it is possible to calculate the radiated field from 
the sphere-surface system in a particular direction given 
by the wave vector q', 
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If we consider a direction other than the specular di- 
rection, the radiated field comprises two contributions: 
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where E ra< j, sp here represents the field radiated directly by 
the sphere, and E ra d jSU rfacc corresponds to waves having 
a final scattering event from the surface. 



III. RESULTS AND DISCUSSION 

In this section, we study the system formed by a gold 
sphere and a gold surface by means of the theoretical 
method described above. The dielectric function of gold 
is described using the tables reported in Ref. [32]. We 
will consider an incident p-polarized wave with angle of 
incidence Q\ n = 57°. We will mainly concentrate on 
the effects associated with the formation of resonances 
of the tip-surface coupled system. In what follows, we 
will mainly focus on the analysis of the local field en- 
hancement and the scattering cross section for scattered 
p-polarized light. Results will be shown for the obser- 
vation angles O ut = 57° and aw = 90°. It should be 
noticed that similar effects than the ones discussed in the 
present work can be seen for other angles of observation 
and incidence. 



Figure ||(a) shows the local-field enhancement below 
the sphere as a function of energy for different values of 
the tip-sample distance and a sphere radius of 40 nm. 
The peak that is built up as the sphere approaches the 
surface is due to a surface plasmon of the coupled system 
formed by the sphere and the surface.Lj This maximum 
experiences a redshift as the tip approaches the surface. 
The scattered light from the system shows a similar peak, 
as illustrated by Fig. ^(b). 

Figure ||(a) shows the local-field enhancement when a 
larger sphere (radius of 80 nm) is considered. The first 
peak has moved to lower energy (compared to the previ- 
ous case with smaller radius) and a new peak appears at 
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FIG. 2: (a) Near-field enhancement just underneath the 
sphere and (b) differential scattering cross section as a func- 
tion of the energy of the incident wave for a gold sphere of 
radius R=40 nm, and for different values of the distance be- 
tween the tip and the surface, d. The angle of incidence is 
#m = 57° and the angles of observation are 9 ou t = 57° and 
(pout = 90°. 
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higher energy, but the field enhancement does not vary 
substantially over a wide range of energies. The scat- 
tering cross section has a qualitatively different behavior 
[see Fig. ||(b)]. The first peak appears also in the scatter- 
ing cross section. By contrast, for higher energies, there 
is no significant correspondence between the scattering 
cross section and the local field enhancement. 

In order to analyze the different features in 
Figs. H and |[ we study the distribution of the electric 
field near the sphere. Figure |J represents the electric 
field around the 80 nm-radius sphere for the first field 
enhancement maximum in Fig. y(a). This figure, as well 
as other vector plots shown below, displays the electric 
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field just outside the sphere in a vertical plane that also 
contains the incident wave vector q. The fields shown in 
Fig. [| corresponds to the fundamental resonant surface- 
plasmon mode of the tip-surface system. This funda- 
mental mode is also responsible for the peak in Fig. ^(a), 
where a smaller radius is considered. The local field in 
that case shows a very similar pattern. However, for the 
second peak in Fig ||(a), the field pattern is significantly 
different, as can be seen in Fig. The local field shows 
extra features that correspond to a higher mode of the 
tip-surface system. The dipolar character of the funda- 
mental mode makes it fairly efficient in emitting light, 
as shown by the large scattering cross section associated 
with it (see Fig. ^(b)). On the contrary, the higher mode 
has a small scattering cross section associated, which can 
be explained by the more complex distribution of the lo- 
cal field for this mode. 

The formation of such resonances and their energy 
shifts have significant consequences on the "approach 
curves", i.e. when varying the tip-surface distance for 
fixed values of the photon energy. Figure ^ shows the 
local-field enhancement and the scattering cross section 
for a photon energy of 1.9 eV and different values of the 
radius of the sphere. The local-field enhancement in- 
creases when the tip approaches the surface for the range 
of tip-sample distance studied and for all the sphere radii 
under study [Fig. 0(c)]. (We only show results for tip- 
surface distances larger than 0.3 nm, since for very short 
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FIG. 3: (a) Near-field enhancement just underneath the 
sphere and (b) differential scattering cross section as a func- 
tion of the energy of the incident wave for a gold sphere of 
radius R=80 nm, and for different values of the distance be- 
tween the tip and the surface, d. The angle of incidence is 
#i n = 57° and the angles of observation are ou t = 57° and 
(pout = 90°. 



FIG. 4: Real part (a) and imaginary part (b) of the electric 
field on the sphere surface for R — 80 nm, E = 1.32 eV 
and d = 1 nm, which corresponds to the first maximum of 
the near- field enhancement in Fig. [](a). The displayed fields 
have been calculated just outside the surface of the sphere in 
a vertical plane that contains the incident wave vector q. The 
absolute magnitude of the depicted fields is about 325 times 
stronger than the incident field as can be seen in Fig. fa. 



5 



distances the description of the system withj-a macro- 
scopic dielectric function is no longer justified £3) 

Figure ^|(a) illustrates the behavior of the differential 
scattering cross section as a function of tip-surface dis- 
tance. It presents a maximum. It is worth noticing that 
this maximum does not imply a maximum of the lo- 
cal field enhancement, which keeps increasing when the 
sphere approaches the surface. The corresponding fea- 
ture in the local field enhancement is a shoulder, rather 
than a maximum, as can be seen in Fig ^|(c) for sphere 
radius R=40 nm. The tip-surface distance at which the 
maximum of the scattered field appears increases with 
the radius of the sphere and the peak becomes wider. In 
addition, some extra structure appears at very short dis- 
tances [see fig. |](b)], especially for the larger radii. We 
will analyze the physical origin of the different features 
by looking in detail at the local field distribution. 

Figure [7] displays the local field around the sphere 
for the maximum of the scattering cross section for a 
sphere with radius 40 nm. This figure shows that the 
peak corresponds to the fundamental surface plasmon 
mode of the tip-surface system. The approach curves, 
for small spheres were previously calculated by Girard,E-j 
who modeled the sphere by an electric dipole and ob- 
tained peaks that can be related to the fundamental sur- 
face plasmon mode. For the larger radii under study, 
the fundamental mode resonance is already shifted to 
energies lower than 1.9 eV, the photon energy considered 
in Fig. ^. The behavior of the scattering cross section 
as a function of tip-surface distance is in this case con- 
trolled by the tail of the resonance peaks shown in Fig. 



o 

-0.2 
-0.4 
-0.6 
-0.8 
-1 
-1.2 



en 



-1.4 


-0.2 
-0.4 
-0.6 
-0.8 
-1 
-1.2 
-1.4 



- (a) \ 

; '""'"'n 




t i 
- (b) : 

; "'■■■■..,„/ 


\ : 



-1.5 -1 



-0.5 0.5 1 

X/R 



1.5 



FIG. 5: Real part (a) and imaginary part (b) of the electric 
field on the sphere surface for R = 80 nm, E = 1.9 eV and 
d = 1 nm, which corresponds to the second maximum (with 
a magnitude ~ 300 times that of the incident field) of the 
near-field enhancement in Fig. ^(a). 
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FIG. 6: Differential scattering cross section (a) and (b), and 
near-field enhancement just underneath the sphere (c), as a 
function of the tip-surface distance d, for a photon energy of 
1.9 eV and sphere radius R of 40 nm, 60 nm, and 80 nm. The 
angle of incidence is t9j n = 57° and the angles of observation 
are t9 out = 57° and <j> out = 90° . 



|3|. This explains why the peaks in the scattering cross 
section as a function of tip-sample distance displayed in 
Fig. |](a) and calculated for larger sphere radii (60 nm 
and 80 nm) are much wider than the one calculated for 
the 40-nm sphere. At a much smaller tip-surface dis- 
tance, i.e. about 0.5 nm, a narrow peak is present for the 
i?=80 nm sphere [Fig. ^(b)]. The local field distribution 
displayed in Fig. with its characteristic node, shows 
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that this peak is associated with the higher mode reso- 
nance of the tip-surface system. The tip-sample distance 
at which this higher-mode resonance appears increases 
with increasing sphere radius and photon energy. 

When studying the light scattered by an isolated 
sphere one finds that for the larger sphere (R=80 nm) 
there is a Mie resonance at energies near 1.9 eV. When 
considering the sphere-surface system, the Mie resonance 
of the isolated sphere implies that, for the 80 nm-radius 
sphere, the amount of s-polarized scattered light be- 
comes larger than the p-polarized light. Nevertheless, 
the amount of s-polarized scattered light is almost inde- 
pendent of the tip-sample distance, except for distances 
shorter than 1 nm, where some features appear that could 
be related to the formation of the higher resonant mode 
discussed above. 
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FIG. 7: Real part (a) and imaginary part (b) of the electric 
field on the sphere surface for R — 40 nm, E = 1.9 eV and 
d — 1.23 nm, which corresponds to a peak of scattering cross 
section in Fig. ^(b) . The local field is about 400 times stronger 
than the incident field in this case as can be seen in Fig. H(c) . 



It is also interesting to analyze the approach curve 
for photon energies with no significant resonance effects. 
Figure ||(a) shows the local-field enhancement for a pho- 
ton energy of 2.4 eV. In this "off-resonance" situa- 
tion, the local-field enhancement presents a maximum 
as a function of the tip-sample distance. This maximum 
appears between d=l nm and d=3 nm for the sphere 
radii studied. The maximum of the local field enhance- 
ment manifests itself also in the scattering cross section 
through a very small peak at similar tip-surface distance 
[see Fig. ||(b)]. However, the scattering cross-section 
maximum appears at much larger tip-surface distances, 
between d=10 nm and d=20 nm. 

In order to gain physical insight about the origin 
of the different maxima in the off-resonance situation 
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FIG. 8: Real part (a) and imaginary part (b) of the electric 
field on the sphere surface for R — 80 nm, E = 1.9 eV and 
d — 0.48 nm, which corresponds to a peak of scattering cross 
section in Fig. |(b). From Fig. |(c) we see that the local field 
in this case is about 1000 times stronger than the incident 
field. 



[E=2A eV), we study the convergence with the number 
of multipoles taken into account in the calculation. Wc 
analyze the case of a sphere radius of 40 nm, for which 
the maxima of enhancement and scattered amplitude ap- 
pear at well separated values of tip-surface distance. The 
local-field enhancement requires a high number of multi- 
poles to converge, in this particular case it was necessary 
to consider multipoles of order I — 20 to arrive at a sat- 
isfactory convergence. Differently, convergence for the 
scattering cross section is achieved with a low number of 
multipoles. Retaining multipoles up to I — 2, i.e. dipole 
and quadrupole modes of the sphere, we obtain most of 
the features of the scattering cross section, in particular, 
its main maximum. This fast convergence implies that a 
simple model for the sphere-surface system can be ade- 
quate for the estimation of the scattering off the system 
when no resonant behavior is involved. In particular, 
the often used model consisting of a dipole represent- 
ing the sphere and its image dipole that together have 
a quadrupole character might be suitable in this situa- 
tion. If we are interested in the local-field enhancement, 
even in this off resonance situation, we need to consider 
a more sophisticated model including higher multipoles 
of the sphere. 

On the other hand, for a photon energy of 1.9 eV (the 
"on-resonance" situation for which cavity modes can be 
built up), it was necessary to take into account a high 
number of multipoles in order to converge both the scat- 
tering cross section and the local field enhancement. For 
example, for a sphere radius of 60 nm, for reproducing 
the main wide peak of the scattering cross section it was 
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FIG. 9: (a) Near-field enhancement just underneath the 
sphere, and (b) differential scattering cross section as a func- 
tion of the tip-surface distance d, for a photon energy of 2.4 eV 
and sphere radii of 40 nm, 60 nm, and 80 nm, respectively. 
The angle of incidence is 9 ln — 57° and the angles of observa- 
tion are 9 out = 57° and </w = 90° . 



enough to retain multipoles up to around I — 10. How- 
ever, we needed to take into account multipoles up to 
the order of I — 70 to arrive at convergence for the lo- 
cal field enhancement and to reproduce the peak of the 
scattering cross section at short sphere-surface distance 
associated with a higher resonant mode. For the range 
of energies where tip-surface plasmon modes can be built 
up, a simple model is inadequate to describe the system. 

Cavity or gap modes h aye been reported in the context 
of light emission in STMESlj and in samples consisting of 
metallic nanoparticles placed very close to a metallic sur- 
face.!^ In both cases, the distance between the metallic 
tip or nanoparticle and the metallic surface can be shorter 
than 1 nanometer. For SNOM, little has been done in 



connection with gap modes. In the context of aperture- 
less SNOM, the approach curves have bee 
tally studied in detail by several groups.! 
A usual effect is the appearance of a peak or maximum 
of the signal when tip approaches the sample, especially 
when both the tip and the sample are made of the same 
metaLEj This peak, which is sensitive to the polarization 
and the dielectric functions of the tip and the sample, 
can be interpreted as a manifestation of the fundamental 
cavity mode discussed in this article. In a recent arti- 
clejid Andre et al. have reported experimental results on 
apertureless SNOM with a STM that might be related to 
the formation of cavity modes between the tip and the 
surface. Particularly, different behaviors of the signal as 
a function of the tip-substrate distance have been mea- 
sured, with the appearance of minima and maxima as the 
tip is removed from the substrate. There are some qual- 
itative similarities with the results shown in Fig. plsup- 
porting the claim that the effects observed in Ref.^l] are 
associated with cavity modes. Nevertheless, the experi- 
mental setup in Ref. [ll] is different than the one consid- 
ered in the present work, therefore the comparison should 
be established with caution. 

So far, we have considered sphere-surface distances 
much shorter than the wavelength. For sphere-surface 
separations on the order or larger than the wavelength of 
the incident wave, both the scattering cross-section and 
the near-field enhancement oscillate as a function of the 
tip-surface distance. These undulations are due to the 
interference of the incident wave and the wave reflected 
by the surface. At this range of sphere-surface distance, 
there is no significant coupling between the sphere and 
the sample. 



IV. CONCLUSIONS 

By means of a detailed theoretical framework, which 
includes retardation, higher multipoles of the tip, and 
multiple scattering between the tip and the sample, we 
have studied the effects associated with the tip-sample 
coupling in apertureless SNOM. In particular, we have 
analyzed the implications of energy shifts of the tip- 
induced surface plasmons or cavity modes as the tip ap- 
proaches the surface. There is no exact correspondence 
between the local field enhancement and the scattering 
cross section. As a function of the tip-surface distance, 
the scattering cross section present maxima related to 
the cavity modes, while the local field enhancement keeps 
increasing as the tip approaches the surface. For a pho- 
ton energy for which resonant effects are negligible, both 
the local field enhancement and the scattering cross sec- 
tion present maxima but at well separated values of the 
tip-surface distance, in contrast to what is sometimes as- 
sumed. In some cases, the maximum of the scattering 
cross section appears at tip-surface distances of 15 nm 
while the maximum local-field enhancement is at a tip- 
surface distance of approximately 2 nm. 
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We have also studied the convergence of the results in 
terms of the number of multipolcs included in the cal- 
culation. In general, one must retain a high number of 
multipolcs in order to arrive at convergence. Only for 
energies for which no cavity modes are built up, a sim- 
ple model (retaining dipole and quadrupole modes of the 
sphere) for the sphere-surface system can be appropri- 
ate for estimating the light scattered off the system, but 
insufficient for the local held enhancement. When tip- 
surface plasmon modes are present, the simple model is 
inadequate and a high number of multipoles should be 
taken into account. 

Finally, we would like to stress that the present results 
are-also significant for surface-enhanced Raman scatter- 
ing^ and the emerging area of tip-enhanced Raman spec- 



troscopy.E2lc-SE3 The formation of cavity modes between 
the tip and the sample and the associated field enhance- 
ment can play a significant role.cB 
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